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We present a determination of the strong coupling constant using ATLAS inclusive jet cross section data at 
-\fs = 7TeV, with their full information on the bin-to-bin correlations. Several procedures for combining the sta- 
tistical information from the data inputs are compared. The theoretical prediction is obtained using NLO QCD with 
non-perturbative corrections. Our determination uses inputs with transverse momenta between 45 and 600 GeV, the 
running of the strong coupling being also tested in this range. Good agreement is observed when comparing our result 
with the world average at the Z-boson scale, as well as with the most recent results from the Tevatron. 
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1. Introduction 

In the Standard Model of particle physics (SM), the 
strong coupling constant a s is one of the fundamental 
parameters. Testing the energy dependence (running) of 
a s over a wide range provides an implicit test of QCD 
and probes potential effects from "New Physics". 

We present a determination of a s using the ATLAS 
inclusive jet cross section data [1], allowing for the first 
test of the running of a s up to the TeV scale. For a 
given a s value, the perturbative NLO QCD prediction 
of the inclusive jet cross section is computed, and is 
then corrected for non-perturbative effects. The pertur- 
bative computation uses the NLO QCD proton parton 
density functions (PDFs), as obtained by various groups 
from other independent data sets. We first determine 
a s (M|) using the measured cross section in each data 
bin, and then propagate the experimental uncertainties 
and bin-to-bin correlations. Finally, an unbiased aver- 
aging procedure is used and the uncertainties are propa- 
gated to the corresponding result. 

The results presented here are based on the studies 
described into detail in |2||. 



2. Input Data 

In this study we exploit the unfolded double- 
differential ATLAS inclusive jet cross section data, in 
bins of the absolute rapidity (|y|) and as a function of 



the jet transverse momentum (pi), with the statistical 
and systematic uncertainties, and their correlations |Q]]. 

This measurement, performed for anti-A;, jets [3] with 
a distance parameter R=0.4 and R=0.6 respectively, 
covers a large phase-space region going up to 4.4 in 
rapidity and from 20 GeV to more than 1 TeV in pj. 
The detailed analysis of the uncertainties and their cor- 
relations is one of the important achievements of this 
measurement. 

The dominant uncertainty in this measurement is due 
to the jet energy scale (JES) calibration. Other smaller 
systematic uncertainties are due to the luminosity mea- 
surement, multiple proton-proton interactions, trigger 
efficiency, jet reconstruction and identification, jet en- 
ergy resolution and deconvolution of detector effects. 
The amplitudes of these (asymmetric) uncertainties and 
their (anti-) correlations between the various |y| and 
pi bins are provided using a set of 87 independent nui- 
sance parameters. Depending on the rapidity bin, the 
systematic uncertainties account for 20 - 60% at low 
pi, and 20 - 40% at high pi, while they are smaller for 
intermediate jet transverse momentum values. 

The statistical uncertainties and their (anti-) correla- 
tions are provided through a covariance matrix for each 
rapidity bin. They are in general much smaller than the 
systematics, except for the high pi region, where they 
can be larger. 

It has been pointed out in that, when performing 
the comparison between data and the theoretical pre- 
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diction, somewhat larger differences are observed for 
R=0.4 than for R=0.6 jets. In our study, we do not per- 
form a simultaneous treatment of the two jet sizes, be- 
cause the correlations of the uncertainties of these mea- 
surements (expected to be large) are not available. We 
use the measurement of jets with R=0.6 for determining 
our nominal result, while the comparison with the result 
using jets with R=0.4 provides a systematic uncertainty. 

3. Theoretical Prediction 

Theoretical predictions for the inclusive jet cross sec- 
tions are calculated in perturbative QCD at next-to- 
leading order (NLO) using the NLOJET++ 0] pro- 
gram. The effects of hadronisation and underlying 
event activity on jet production have been studied using 
Monte Carlo generators. Each bin of the perturbative 
cross section is multiplied by the non-perturbative cor- 
rections to obtain the prediction at the particle level that 
can be compared to the data. 

To allow for fast and flexible evaluation of the 
cross section using different values of a s (exploiting 
the a s scans performed in the PDF fits), as well as 
for evaluation of the PDF and scale uncertainties, the 
APPLGRID H software was interfaced with NLO- 
JET++ in order to calculate the perturbative coefficients 
once and store them in a look-up table. The nominal 
theoretical prediction is derived with the CT10 [6] NLO 
parton density functions, the corresponding uncertain- 
ties (eigenvectors) being propagated too. Uncertain- 
ties due to the choise of the PDF set are evaluated by 
comparing the nominal results with the ones obtained 
with MSTW 2008 H, NNPDF 2.1 flU and HERA- 
PDF 1.5 fill- 
in our analysis we have used the non-perturbative 
corrections and uncertainties evaluated in |l|], using 
leading-logarithmic parton shower generators with their 
various tunes. For R = 0.4 the correction factors are 
dominated by the effect of hadronisation and are ap- 
proximately 0.95 at jet pi = 20 GeV, increasing closer 
to unity at higher pj. For R = 0.6 the correction factors 
are dominated by the UE and are approximately 1 .6 at 
jet pj = 20 GeV decreasing to between 1.0-1.1 for jets 
above p T = 100 GeV. 

4. a s Evaluation Procedure 

The determination of a s (M 2 ) goes through exploiting 
the experimental information in individual (pj; \y\) bins 
of the jet cross section, together with the propagation 
of the (asymmetric) uncertainties and bin-to-bin corre- 
lations. 



An a s value is obtained in each bin of the measure- 
ment, using the theoretical prediction and the experi- 
mental cross section. The theoretical prediction pro- 
vides, in each (pj; \y\) bin, a one-to-one mapping be- 
tween a s and the inclusive jet cross section. The nom- 
inal value of the measured cross section is mapped by 
this procedure to the nominal a s value in the corre- 
sponding bin. 

A series of pseudo-experiments, each of them in- 
cluding a complete set of cross section values in all 
the (pj; \y\) bins, are generated using the experimen- 
tal input described in section [2] The statistical fluc- 
tuations are generated using the information provided 
by the covariance matrices. Each nuisance parameter 
for a given component of the systematic uncertainties 
is treated as 100% correlated between all the bins. The 
various nuisance parameters are treated as independent 
with respect to each other. The distributions of the ex- 
perimental uncertainties are modelled using asymmetric 
Gaussian shapes, centered at zero, with 50% probabil- 
ity of generating pseudo-measurements at negative or 
positive values. Other possible models of these distri- 
butions have been tested, yielding similar results. For 
each pseudo-experiment, the same procedure as for the 
nominal measurements is applied to obtain an a s value 
in each bin. In view of the study of various possible av- 
eraging procedures, it is interesting to determine (after 
a small symmetrisation of uncertainties) the covariance 
matrix of the a s values. 

In order to avoid potential biases introduced by some 
averaging procedures, several options for performing 
the combination of the inputs in the various (pj; \y\) bins 
have been tested. All of them aim at obtaining an aver- 
age value of the strong coupling constant (a^ v ). 

The first procedure consists of a simple average, 
where all the input a s values have the same weight. 
The 2" d procedure which has been tested consists of a 
weighted average, where the weights of the input a s val- 
ues are proportional to the inverse of their squared total 
uncertainties. Although the final uncertainties on these 
two averages might not be optimal, the weights of the 
input a s values are well behaved (i.e. they are all in the 
[0; 1] interval, with the sum equal to unity). The 2 nd 
method also has the advantage that more precise contri- 
butions get larger weights. 

The 3 rd averaging procedure consists of the minimi- 
sation of a "standard"^ 2 with correlations 

X 1 = (a^ om - of) ■ CT 1 ■ (o£ om - off, (1) 

with respect to a^ v . Here, af om is the vector of nominal 
a s values (obtained in individual (pj; \y\) bins) entering 
the average, C is their covariance matrix, while a^ v is a 



2 



vector containing values equal to af* (to be fitted) and 
having the same size as a^ om . 

It has been shown in [11] that this approach, al- 
though very commonly used, can yield biased average 
values. Actually, in presence of (not very well under- 
stood) strong correlations among the inputs, the average 
value can be even outside the range of the input values, 
which seems unacceptable. Indeed, the minimisation of 
the function in Eq. [TJ with respect to a^ v (equivalent 
to the minimisation of the variance of the output of a 
weighted average, cf. Gauss-Markov theorem), yields 
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where 1 is a vector with all the entries equal to unity 
and the same size as a^° m . Just as in the two previous 
methods, the average resulting from the^ 2 minimisation 
is a weighted mean of the input values. However, in 
presence of strong correlations, these weights can even 
be negative or larger than unity. 

The explanation of the bias in this procedure, as pro- 
vided in jllll . makes responsible the input co variance 
matrix because of "the linearisation on which the usual 
error propagation relies". In 11211 . the bias was ex- 
plained by the use of the ill-behaved inversion of the 
covariance matrix, as well as the treatment of all the 
uncertainties as absolute (i.e. corresponding to additive 
effects), whereas at least some of them should be treated 
as relative uncertainties (i.e. corresponding to multi- 
plicative effects). In yO we have pointed out one addi- 
tional problem: this combination procedure fully relies 
on the knowledge of the bin-to-bin correlations (difficult 
to check in a particular measurement), in order to derive 
"optimal weights", providing the smallest uncertainty 
on the average. This yields too optimistic results, as 
the uncertainties on the shape and correlations for each 
component of the systematic uncertainties should also 
be taken into account. 

The three combination procedures described above, 
make use of weighted averages. They distinguish be- 
tween each other by the way how the weights of the 
inputs in the average are obtained. Another possibility 
that we have considered for performing the combination 
consists in a geometrical average. 

5. Experimental results 

As explained in the previous section, a nominal 
a s value has been determined in each (pj; \y\) bin, 
while the statistical and systematic uncertainties, to- 
gether with their correlations, have been propagated us- 
ing pseudo-experiments. For very low and very large 



pj values, the total data uncertainties are relatively 
large, and a fluctuation of the cross section by ±l<x (or 
more) yields a s values for which the theoretical predic- 
tion is not reliable. It is for this reason that we discard 
the pi bins with large experimental uncertainties (hence 
small weights) in the final a s average computation, and 
we estimate a systematic uncertainty due to this choice. 
A number of 42 (pj\ \y\) bins are used in the computa- 
tion of the various averages. 

The precision of the weighted average (WA) com- 
puted in the individual \y\ intervals of the experimental 
measurement |l[], is similar to the precision of the val- 
ues obtained in the corresponding (input) pj bins |01- 
This is due to the strong correlations between (and the 
similar size of) the uncertainties in the pj bins used for 
the average inside a given \y\ bin. Performing the com- 
bination of all the 42 {pj; \y\) bins, one obtains 

ff r =o.ii5i™7. (3) 

The experimental uncertainty of this average is about 
30% smaller compared to the one obtained in the most 
precise individual \y\ bin (i.e. the central bin). This gain 
in precision is due to the rather small correlations of 
the systematic uncertainties, between the central and the 
forward region. The results (central values and uncer- 
tainties) obtained using the weighted average, the sim- 
ple average and the geometrical mean are consistent be- 
tween each other, as expected for a set of inputs with 
similar uncertainties. 

The result of the weighted average, together with the 
covariance matrix of the uncertainties for the a s val- 
ues, have been inserted into Eq.[JJ yielding a^- 2 /dof = 
0.54 (for 41 degrees of freedom). This represents an 
implicit test of the RGE, used for the perturbative pre- 
diction as well as the fits for determining the PDFs used 
here. Very reasonable values for the^ 2 /dof are obtained 
in the individual \y\ bins too. Taking into account the 
bin-to-bin correlations is very important for the deter- 
mination of these x 1 values. Ignoring the correlations 
would actually yield x 1 values smaller by more than one 
order of magnitude. 

The procedure based on the minimisation of a^- 2 with 
correlation (see Eq. [TJ yields 



or,™" =0.1165 



+0.0033 
-0.0033- 



(4) 



Compared to the previous averaging procedures, this re- 
sult has a smaller experimental uncertainty and slightly 
smaller ^f 2 /dof (0.53), as expected. However, the 
central value obtained here is questionable, since the 
weights of this average strongly rely on the exact knowl- 
edge of the bin-to-bin correlations and are not well be- 
haved. Indeed, a large fraction (almost half) of the 
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weights of the individual pj bins in this average are 
smaller than zero. Looking at the similar averages com- 
puted in individual |y| bins, they are in the range between 
0.1100 and 0.1243, several of them being outside the 
range of the corresponding input values. This is due to 
the same bad behaviour of the weights, some of which 
are negative, while some others are larger than one. 

It is due to the issues observed for the minimisa- 
tion of a x 2 with correlation, compared to the well be- 
haved weighted average, that we choose to use the value 
obtained with the weights given by the inverse of the 
squared total uncertainties, for our nominal result. 
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Figure 1: Weighted a s average (green band) evolved to the corre- 
sponding pi scale, together with the values obtained in all the (pj; 
\y\) bins used in the combination. The blue vertical line indicates the 
highest pj value used in the Tevatron a s (M|) determination, from the 
inclusive jet cross section. 

The nominal result for a s evolved by N 3 LO RGE run- 
ning, its uncertainty, and the comparison with the in- 
puts used in the combination are shown in Fig. Q] We 
also indicate the highest pj value (145 GeV) used in the 
a s (M^) determination, from the inclusive jet cross sec- 




tion, at Tevatron 111 311 . Our combination procedure uses 
a pi range going up to four times higher. We also use 
inputs from the forward region, up to \y\ = 4.4, while 
the Tevatron determination is restricted to |y| < 1 .6. 

Results with a similar precision are obtained using 
anti-A:, jets with R=0.4. The central value of the corre- 
sponding weighted average is, however, shifted down- 
wards with respect to our nominal result (from anti- 
k, jets with R=0.6), by about 0.0060. This is expected, 
due to a systematic difference seen between the com- 
parisons of the two (strongly correlated) measurements 
with the corresponding theoretical predictions (see sec- 
tion The difference between the a s determinations 
with the two jet sizes is actually our largest uncertainty 
and an improvement in the understanding of these dif- 
ferences is desirable. 



6. Conclusions 

We have performed a determination of the strong cou- 
pling constant at the Z scale, using the ATLAS inclusive 
jet cross section data. Our final result accounts for J^] : 

a s {M 2 z ) = 0.1151 ± 0.0001 (stat.) ± 0.0047 (exp. syst.) 
±0.0014 (p T range) ± 0.0060 (jet size) 
+0 00 lt (scale)«2 (PDF choice) 



-0.0011 



+0.0010 (PDF eig.)!o:oo34 < NP corrections), 

where all the statistical and systematic uncertainties are 
included. Our value is in good agreement with the 
latest (preliminary) update of the a s (M^) world aver- 
age (0.1 183 ± 0.0010) lQ, as well as with the latest re- 



sult from a hadron-hadron collider (0.1161+°^) |13H . 



Although our result is less precise, it includes for the 
first time the measurements of the inclusive jet cross 
section up to 600 GeV. The running of a s has also been 
tested, in the pj range between 45 and 600 GeV, and 
no evidence of a deviation from the QCD prediction has 
been observed. 
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